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Abstract 


This paper describes the application of neural network adaptive wavelets for fault diagnosis of space station power system. 
The method combines wavelet transform with neural network by incorporating daughter wavelets into weights. Therefore, the 
wavelet transform and neural network training procedure become one stage, which avoids the complex computation of wavelet 
parameters and makf* the procedure more straightforward. The simulation results show that the proposed method is very 
efficient for the identification of fault locations. 

INTRODUCTION 


Wavelet transform (WT) shows promise for both signal (or image) representation and classification. It has been successfully 
used in many areas <a»ch as optical computation, seismics, acoustics, and mechanical vibrations (Szu 1992, Rioul 1991, and 
Chui 1992). The representation and classification both can be viewed as feature extraction problems in which the goal is to 
find a set of daughter wavelets that either best represent the signal or best separate various signal classes in the resulting 
feature space. For these applications, wavelet analysis overcomes the limitations of Fourier method by employing analysis 
functions that are local both in time and frequency and allows the convenient representation and c l ass ifi cation of duration 
signals. Recently, several papers have been presented for power system transient signal analysis (Riberio 1994, Santoso 1995, 
Wilkinson 19%, Pillay 19%, and Momoh 19%). These works focused on the representation of transient signal. An 
optimization procedure is needed for the calculations of wavelet transform coefficients (3j, bi). This procedure may be time 
consuming. On the other hand, the best set of wavelets for representation will not necessarily be the same as the best set for 
classification. As a new tool, wavelets have not been applied to the classification problem of power system. 

More recently, Artificial Neural Network (ANN) techniques have been used to solve fault diagnosis problem in power 
system. The use of neural network based scheme for fault detection and classification is only in its first stage (Chowdhury 
19%, and Momoh 19%). In these works, a pre-processing of fault signals by Fast Fourier Transform (FFT) is used It is a 
two stage process. The first stage is feature extraction, the second stage is the training of ANN. The disadvantage of this 
scheme are: (1) The procedure may be time consuming, and (2) FFT technique can be applied with success to steady state 
phenomena but short-time events require different mathematical tools to aid the analysis. 

In this paper, we discuss how wavelet transform coefficients can be adaptively computed in the training stage of ANN. Each 
daughter wavelet is represented by a weight from input layer to first hidden layer. The neurons of the first hidden layer differ 
from the neuron of other hidden layers. The inspiration functions of neurons in first hidden layer are line a r functions, so that 
the first hidden layer introduced for wavelet transform can be compressed to the second layer. The wavelet transform 
coefficients are automatically computed in the training stage of neural network. This neural network adaptive wavelet scheme 
is applied to the fault diagnosis of space station power system. 

The space station power system will provide electrical power to the future space station by converting the solar energy 
collected from the sun into electrical energy. An AC/DC system architecture is recommended for this special power system. 



Figure 1 shows the system structure. The system includes a generator, AC network, an AC/DC converter, DC network, 
DDCUs and loads AC system is a generation system which transfers solar energy to electrical energy, while DC system is a 
distribution system which provides electrical power to all kinds of loads. The generation and distribution systems are 
connected by an AC/DC converter. Since the time constant of the distribution system is very short, the transient process of the 
distribution system is very short. Wavelet is, therefore, an efficient choice. The fault diagnosis is carried out for distribution 
system by the proposed neural network adaptive wavelets. The output current of DDCU is recorded for fault diagnosis. 

The simulation results show that the proposed method is very efficient for the identification of the distribution system faults 
at different locations. The misclassification band for most faults is about 10% -15% of line length from the sending end bus. 



FIGURE 1: Space S tation Po wer Sy stem 

NEURAL NETWORK ADAPATIVE WAVELET MODEL 

Different models are used for the representation and classification of signals. A short discussion on representation model and 
classification model is as follows. 


Representation Model 


A signal s(t) can be approximated by daughters of a mother wavelet h(t) according to 
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where the {w k }» {bjjand {a*} are weight coefficients, shifts, and dilation for each daughter wavelet. Generally, an 
optimization procedure is used for the calculations of these parameters, which may be time-consuming. In order to speed up 
the procedure, the approximation of equation {!) can be expressed as the neural network of Figure 2, which contains wavelet 

nonlinearities (f i = h (- — — )} in the artificial neurons in hidden layer rather than the standard sigmoidal nonlinearities. The 

network parameters fw k }, {t\} and {a*} can be optimized by minimizing an energy function. This optimization procedure is 
the same as the procedure for the training of ANN. Therefore, the two optimization procedures can be combined together, i.e., 
the optimization procedure fee wavelet transform is integrated into the optimization procedure of ANN training. The least - 
mean - squares energy function is employed for signal representation. 
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FIGURE 2. Representation Neural Network 



FIGURE 3. Classification Neural Network 


Classification Model 


The extraction of signal features is the vector inner products of a set of wavelets with 
then be input to a classifier. The combined classifier and wavelet feature detector is given 
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the input signal, 
by 

(3) 


These features can 


where v B is the output for the n 01 training vector Sn(t) and a(z) = ^— r- This classifier can be depicted as the neural 
network of Figure 3, which uses wavelet weights rather than the wavelet nonlinerities of as in Figure 3. The lower part of 
Figure 3 produces inner products of the signal which has form of 'Ls(i) h( — ~) and wavelets with the first wavelet on the 

left and the nk* wavelet on the right 


Figure 3 shows two layers of weights. /„ = h(-^-) represents the wavelet transform coefficient, while the classification 

^ J 

parameters {w k }, {a^Jand {bt} can be optimized by minimizing the error between the calculated and desired outputs. 
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where d B is Uie desired classifier output for s„(t). It can be set such dial do = 1 for one class and zero for the other. The 
classification neural network can be simply represented by Figure 4 (a). The inputs of the ANN are the classified signal values 
at time intervals (t=l, 2, .... T). A network with more output elements is shown in Figure 4 (b). The input signals of ANN are 
the same as the single case. The desired classifier outputs (dj, d 2 , .... d„) are set such that 


( d, = 1, if case i is I me 
\ d t = 0, if case i is not true 
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FIGURE 4. Simplified Classification Neural Network 


METHODOLOGY 


From the models developed above, it can be seen that only one stage is used for ANN training and extraction of signal 
features. The feature extraction is integrated in the traiuing process. The algorithm is the same as the standard ANN training 
algorithm. In our study, the back propagation algorithm (El-Sharkawi 1996) is used. 


According to the back propagation algorithm, the network’s weights {w s } are adjusted to minimize the error between desired 
output and the calculated output. The gradient descent algorithm adapts the weights according to the gradient error, i.e.. 


dE dV. dE 

9 dv dw dw 


(5) 


Specifically, the error signal is defined as 



Equation (5) can be rewritten as 

A W'tsfiSjV, (6) 

Where p is an adaptability factor and ji > 0. 


The algorithm for the training of the adaptive wavelets and ANN is shown in Figure 5. The algorithm consists of the 
following steps. 


(1) Initialization; 

(2) Input classified signal S(l), S(2), .... S(t), and the desired classifier outputs d( (i=l, 2, ...n); 

(3) Calculate the outputs based on the current weights; 

(4) Calculate the error between calculated output and desired output; 

(5) Error check; If convergence is indicated, go to step (6). if not, calculate gradient error and modify weights, go to step (3); 

(6) Output weights and form ANN for fault diagnosis; 

When the ANN is trained, the on line faulted signal can be inputted into the ANN, the output of ANN is the fault diagnosis 
result. 
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FIGURE 4. Flow Chart of ANN Adaptive Wavelet Fault Diagnosis 
THE APPLICATION IN FAULT ANALYSIS 
Test system training data and testing data 

The space station power system architecture is shown in Figure 1. The system data is given in Table 1. System fault locations 
are classified into eleven patterns, which are listed in Table 2. For each case, fault is located along the line, specifically 25%, 
50%, 75%, 100% from the sending bus to receiving bus. A total number of 44 cases were generated for the training of neural 
network. Another 44 faulted cases located at the 10%, 30%, 60%, and 80% of line length are generated for the test of ANN. 

TABLE 1. Space Power System Data 

DDCU Parameters 

K»=200.0, V*=160.0(v), V out =120.0(v), 11=95% 

Input Circuit: 

La=0.704(G), L u =0.5(£2), C„=0.25(Q) Ca=0.2(£2), R„=0.001(G), Ro=0.001(0) Ro=0.005(G). R W =0.005(Q) 

Output Circuit: 

1^=0.005(12), C 0 , =0.00 15(G), C o2 =0.00 15(G) R oi =0.5(n>, R^O.OKG), Rq 3=0.3(G) R^gQ.OOKfl) 

Load data (W): 

P«=250.0, P k =500.0, P, 3=750.0, P M = 1000.0 

P 8 3=1600.0, P 86 =750.0, P 87 =500.0, P, 8 =250.0 

Line data (fl): 

R*=0.001, L4«=0.0002, R$ 7 ,=0.0376, L«7i=0.00136, R« 72 =0.0376, L« 72 =0.00136 

R 71g , =0.00416, L7,8s=0.000051, R 7 , g2 =0.00378, L 7 „ 2 =0.000034, R 71g3 =0.00278, 1^,83=0.000034, R 7 ,84=0.00178. L,, 84=0.000034 
R 72 , 5=0.001 16, L 72 85=0.000051, R 72g6 =0.00278, L 7286 =0.000034, R 72 ; 7 =O.OQ378, 1^87=0-000034, R 72 8,=0.00416, 1^88=0-000034 

System fault analysis 

The fault is assumed to be on a line to ground, and the fault resistance is set to be 0.1 Q. The fault duration is 5ms. The time 
simulation program is used to simulate the system faulted current at the output terminal of DDCU. The DDCU model can be 
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found in (NASA, Lewis Research Center, 1993). Load is modeled as a constant impedance, whereas the line as R-L scries 
circuit model. For each fault, die different locations along the line are simulated to form the training data of ANN. 


TABLE 2. Fault Location Classification 


Fault No. 

Fault along 

Fault No. 

Fault along 

Fault No. 

Fault along 

Fault No. 

Fault along 

F# 1 

line from bus 
4 to bus 6 

F#4 

line from bus 
71 to bus 81 

F# 7 

Line from bus 
7 1 ‘to bus 84 

F# 10 

line from bus 
72 to bus 87 

F#2 ■ 

line from bus 
6 to bus 71 

F# 5 

line from bus 
71 to bus 82 

F# 8 

line from bus 
72 to bus 85 

F# 11 

line from bus 
72 to bus 88 

F#3 

line from bus 
6 to bus 72 

F# 6 

Line from bus 
71 to bus 83 

Fit 9 

line from bus 
72 to bus 86 




Wavelet - ANN Fault Diagnosis 

Weight - based architecture shown in Figure 6 is used for fault diagnosis. The same ANN topology with different weight 
architectures is used to classify different fault pattern. For N patterns of faults, N different architectures are used. Each 
architecture has one output with two desired values of 1 and 0. The weights of architecture (Arch) i (i=l, 2, ..., N) is trained 
in the way that the output of the architecture is close to 1 for F# i and zero for other fault patterns. When the faulted signal of 
F# i is inputted to the fault diagnosis ANN, only the output at architecture i is close to I, the outputs at other architectures are 
close to zero. Therefore, the fault diagnosis is conducted by checking the outputs of the architecture. 


Arch 1 


1: pattern ! fault 
0: no pattern 1 fault 


Faulted 

Current 


Arch 2 


1 : pattern 2 fault 
0: no pattern 2 fault 


Arch i 


,!: pattern i fault 
0: no pattern i fault 



1 : 

0 : 


pattern N fault 
no pattern N fault 


FIGURE 6. Wavelet - ANN Fault Diagnosis Architecture 

Three patterns of faults ( F# 1, F# 2, and F# 4) are used to explain our fault diagnosis scheme. Figure 7 shows the faulted 
currents at DDCU output terminal for faults #1, #3, and #4. The fault is a short-time event, the faulted currents are oscillatory 
and have amplitudes that quickly decay to steady values. The feature meets the conditions of wavelet transform technique. The 
faulted currents at different locations are different. These faulted currents are used for the training of ANN. Figure 8 shows 
the weights for the different faulted currents and the outputs of the corresponding architectures . It can been seen that the 
weights and the outputs of architectures for different faults are quite different which makes it possible for fault diagnosis. 
Since the different fault has different wavelet transform factors, when these transform factors are integrated in the training of 
architecture weights, the different architecture weights will be got for different fault pattern. Different weights will produce 
different outputs, the fault diagnosis can be carried out by checking the outputs of architectures. For a given fault pattern, 
only one output has value close to 1. For example, when faulted current of F# 1 is inputted in the fault diagnosis ANN, the 
output of architecture 1 is close to 1 (1.109088), while the outputs of architectures 2 and 4 are close to zero ( -0.070575 and 
0.06069). Table 3 shows the output values of Wavelet - ANN diagnosis for faulted currents of F#l, F#2, F#3, and F #4. 
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FIGURE 



Table 4 displays llie system fault diagnosis results. Since the system has eleven fault patterns, eleven architectures with 
different weights are used for system fault diagnosis. In this table, True (T) means that the output of ANN according to the 
fault pattern is close to 1 and the fault diagnosis result is correct, while. Failure (F) means that the output of ANN according 
to the fault pattern is far from 1. and the fault diagnosis is wrong. It can be seen that the neural network can identify the faults 
beyond the 10% length of the most lines. When the fault occurs near to the sending end bus, since the resistance of this 
section of line is very small, die faulted current is very close to the faulted current at sending end bus. Therefore, die neural 
network cannot identify the fault correctly. This region is a dead band of the diagnosis. Also, it can be seen that, for line 6-71 
and line 6-72, the fault classifications are not as accurate as the other cases, since the two lines have same parameters, the 
faulted currents are very close when the fault occurs near the sending bus (bus 6). 


TABLE 3. The Output Values of Fault Diagnosis ANNs 



v, 

v 2 

Vj 

KSSSSSI 

F# 1 

1.109088 

-0.070575 

0.07486 

0.06069 

F#2 

0.005727 

0.960932 

0.04939 

0.09996 

F# 3 

0.027072 

0.815518 

0.01037 

0.03019 


0.002880 

-0.11091 

0.11548 

1.15519 


TABLE 4. System Fault Diagnosis Results 


Fault Pattern 

10% 

30% 

60% 

80% 

F# 1 

T 

T 

T 

T 

F# 2 

F 

F 

T 

T 

F# 3 

F 

F 

T 


F# 4 

F 

I 

T 

EBHBEB 

F# 5 

F 

T 

T 

T 

F# 6 

F 

I 

T 

T 

F# 7 

F 

T 

T 


F# 8 

F 

T 

T 

EBBBBB 

F# 9 

F 

T 

T 

EBBBBB 

F# 10 

F 

T 

T 

ebeebbe 

F# 1 1 

F 

T 

T 

EBBBEBB 


* T: Fault diagnosis is correct, F: Fault diagnosis is wrong 
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Weights of Architecture 1 Outputs of Architectures 1, 2, 4 

(a) Fault #1 



Weights of Architecture 2 Outputs of Architectures 1, 2, 4 

(b) Fault #2 



Weights of Architecture 4 Outputs of Architectures 1, 2, 4 

(c) Fault #4 

FIGURE 8. Architecture Weights and Outputs for Different Fault Patterns 
CONCLUSION 

In this paper, we have discussed how wavelet transform coefficients can be adaptively computed in the training stage of 
ANN. This neural network adaptive wavelet scheme is applied to the fault diagnosis of space station power system. The 
architecture weights are trained adaptively for different fault patterns. The test results show that about 85-90% of the line fault 
can be identified by the proposed neural network efficiently. In the proposed scheme, the fault feature extraction stage in 
integrated into the training stage of ANN, which is much helpful for the on-line application. For on-line application, faulted 
signals are directly inputted to ANN for the training instead of a complex processing of all faulted signals. 
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ABSTRACT 


Astro-E is the X-ray satellite to be launched in year 2000 by Inst, of Space & Astronautical Science. 
This report deals with the design and expected performance of the Hard X-ray Detector (HXD), one of 
the 3 experiments aboard Astro-E. The HXD is a combination of GSO/BGO well- type phoswich counters 
and silicon PIN diodes: the two combined will cover a wide energy band of 10 — 700 keV. The detector is 
characterized by its low background of ~ 10 _5 /s/cm 2 /keV and its sensitivity higher than any past missions 
between a few 10 keV and several 100 keV. Combined with the other 2 experiments, a micro- calorimeter 
array (XRS) and 4 CCD arrays (XIS), both with X-ray mirrors, the mission will cover the soft and hard 
X-ray range at a highest sensitivity. 


Keywords: X-ray astronomy, hard X-ray detector, phoswich counter, silicon PIN diode, Astro-E 



1 ASTRO-E PROJECT 


The fifth Japanese X-ray astronomy satellite, ASTRO-E, following Hakucho , Tenma, Ginga, and ASCA , 
is scheduled for launch in year 2000 by the new launcher M-V-4 of Institute of Space and Astronautical 
Science from ISAS Kagoshima Space Center. 1 This satellite will carry three experiments: the Hard X-ray 
Detector (HXD) covering the energy band from 10 keV to 700 keV; a micro-calorimeter array with an X-ray 
mirror (X-ray Spectrometer — XRS); 4 CCDs with 4 X-ray mirrors (X-ray Imaging Spectrometer — XIS). 
The latter two experiments cover the soft X-ray band with the highest energy resolution 2 (micro-calorimeter: 
A E ~ 12eV), or with good energy resolution 3 (CCDs: A E ~ X50eV at 5.9keV), both with imaging capability. 
All 3 experiments combined, Astro- E will become a spectrometer facility covering the energy band from 0.4 
keV to 700 keV with good to moderate spatial resolution (XRS/XRT and XIS/XRT: ~ 1 arcmin., HXD: 
~ 20 arcmin.). 

The total weight and power available for the 3 experiments will be about 800kg and 200 W, respectively. 
The M— V rocket will put the satellite into a near-circular orbit of radius 550km with an inclination of 31 
deg. The 3 scientific instruments, the X-Ray Spectrometer (XRS) with an XRT (focal length ^3. 5m), the 4 
X-ray Imaging Spectrometers (XIS) with an XRT each (focal length ^4. 5m), and the HXD, have adopted 
several new and inovative technologies. The micro-calorimeter array made of 36 HgTe elements will be the 
first such instrument to fly in a satellite and will have the ultimate energy resolution (~12eV) in the soft 
X-ray band. 2 The XRTs will be made of replica foils and will allow 1 arcmin. spatial resolution despite their 
light weight. 4 The X-ray CCDs will have deeper depletion region (~ 80 fim) and their energy resolution 
(fwhm) will reach the theoretical limit for the CCD: below 55eV at Oxygen K-line and 145 eV at 5.9 keV. 3 
The large (~ 20mm x 20mm) 2mm-thick silicon PIN diodes used in the HXD will be the first of this kind and 
will give us an energy resolution about 3keV (fwhm) in the hard X-ray band (~ 10 to 70 keV). The well- type 
phoswich counter itself is a new concept and will use a newly developed fast and high light-yield inorganic 
scintillators GSO(Ce) in its detection part. 5 6 7 These instruments will be prepared by ISAS, Univ. of 
Tokyo, Osaka Univ., Kyoto Univ., Nagoya Univ., Tokyo Metropolitan Univ., RIKEN, NASA/GSFC, Univ. 
of Wisconsin, and MIT. 

The time schedule of the project is as follows: the R/D works have nearly been completed in Japanese 
FY 1995, the Engineering Model (EM) of crucial items with its associated electronic circuitry have been 
produced and tested in the spring of 1996. Based on these tests and R/D works, the design will be frozen 
in the fall of 1996 and production of the Flight Model (FM) will begin early 1997. The FM will then be 
completed in mid- 1998 followed by the final assembly and tests. The launch date is set at present late in 
Japanese FY 1999 or in the winter of calender year 2000. 


2 ASTRO-E HARD X-RAY DETECTOR 


The Astro-E HXD 8 9 has been jointly developed by scientists at Department of Physics, University of 
Tokyo, Institute of Space and Astronautical Sciences (ISAS), Institute of Physical and Chemical Research 
(RIKEN), National Laboratory for High Energy Physics (KEK), and others. It is basically an upgraded 
version of the well-type phoswich counters successfully flown on balloons. 10 11 12 13 14 15 Silicon PIN diodes 
are the important addition introduced to lower the energy coverage as well as to improve the energy resolution 
in the lower energy band. The design and characteristics of the HXD as of June 1996 are described here 
together with its expected performances. 

The HXD detector assembly is schematically shown in Fig. 1. The total weight of this assembly will be 
about 200 kg including the electronic part not shown in the figure. The HXD consists of 16 (=4x4) modular 
units and has an overall photon collecting area of about 330 cm 2 . Each unit is built around a phoswich 
counter made of a fast inorganic scintillator, GSO (Gd 2 Si0 5 :Ce 0.5% mol) 5 6 7 and the BGO active shield 
as shown in Fig. 2. Cosmic hard X-rays are detected as clean-hits if their full energy is deposited in the 
fast scintillator GSO. The active collimation part of the BGO shield forms four deep-wells (320mm deep, 
rv 25mm x 25mm in area), limiting the field of view to ^ 4.3° x 4.3°. The GSO crystals (four per unit, each 
measuring in area 24 mm x 24 mm and in thickness 5 mm) are glued at the bottom surface of the wells. The 
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Figure 1 : The HXD detector assembly: Cross-section (left) and the top view (right). The housing and most 
of the electronic part are not shown here. 
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entire assembly is viewed by a 2 -inch phototube from the exterior surface of the shielding part. 

The silicon PIN diodes, each with an active area of 20 minx 20 mmx 2 mm, are assembled in a unit of 
two layers and placed in the deep BGO well just above the detection part (Fig. 2 ). Two PIN diodes are 
placed per well, making a total of eight per unit. The HXD contains 128 Si diodes in all, achieving a photon 
collecting area of about 230 cm 2 . 

Each group of 4 detector units will be equipped with an independent high voltage supply for the photo- 
tube, and a common DC power supply for 8 PIN diodes. Fig. 3 shows the effective areas of the phoswich 
scintillators and the silicon PIN diodes, in the respective energy ranges of 40 — 700 keV and 10 — 70 keV. 

The 4x4 matrix of phoswich counters are surrounded by 20 units of thick BGO anti-counters for additional 
shielding. Furthuremore, a fine collimator made of phosphor bronze sheet {50/zm thick), is placed inside 
the BGO wells to match the HXD field of view to that of the soft X-ray telescopes (17 x 17 arcmin 2 ). 
This collimeter is expected to reduce the the cosmic diffuse X-ray background that may otherwise become 
a dominant background source for the PIN diodes. In the soft 7 - ray band, background is reduced by 
mutual anti-coincidence among neighboring units. We expect the detector background to be below 1 Q “ 5 
c/s/cm 2 /keV for the Si PIN diodes and ~ 10 -5 c/s/cm 2 /keV for the scintillators. 

The estimated sensitivity will be set by the background due to the radioactivity induced within the 
counters while in the orbit. It will be substantially better if the detector is flown in a lower orbit (eg. that of 
CGRO), where activation is much less, than the ultimate limit the well- type phoswich counter can attain. 
Lowering the background and knowing its origins are crucial in the following respect: We plan to avoid the 
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ON-OFF subtraction by estimating the background accurately (~ a few %) from monitor counts as was 
done in the Ginga LAC observations. 16 

In Fig. 4, the expected 3 a sensitivities to the continuum and to the line emission of Astro-E HXD are 
compared with those of the past and planned satellite missions. The sensitivity of Astro-E HXD for point 
sources will be substantially higher than any other past missions in the energy band between 10 keV to 
several 100 keV. We therefore expect to detect and study many new cosmic hard X-ray point sources. The 
Astro-E HXD, on the other hand, will have only a modest sensitivity to diffuse gamma-ray emissions such 
as the 511 keV line distributed over the Galactic ridge. 


2.1 Well- Type Phoswich Array 

The phoswich counter consists of two kinds of scintillators whose scintillation decay times are distinclty 
different. The faster scintillator is placed in the front as the detection part and the slower one in the back 
as the shielding part. The phototube signals generated purely by the faster scintillator are selected by 
using an appropiate pulse-shape discriminator (PSD). Signals with an appreciable contribution from the 
slower scintillator, eg. those of hard X-rays scattered by the shielding part and those of charged particles 
penetrated through the shielding part, are efficiently rejected (see Sec. 4 . 3 ). This phoswich technique has 
been used for many years in cosmic 7-ray detection. The uniqueness of the “well-type” phoswich counters 
is that the well-shaped shielding part acts also as an active collimator and that each counter filters out, 
quickly and effectively, unwanted hard X-rays down if even only a small energy (50 — lOOkeV) is deposited 
in the collimator and shielding part. 10 12 The detection part (the fast scintillator) being buried deep in the 
active anti- coincidence well also reduces efficiently background due to nuclear activity. The details on the 
development of the well-type phoswich counter are given elsewhere. 10 11 12 13 

Choice of the two scintillation materials becomes important in reducing the background and improving 
energy resolution. For its large effective atomic number and long scintillation decay time, BGO emergied 
as the choice for the shielding part. Through a test on a prototype well, we noted that BGO commercially 
available at that time was contaminated by a radioactive isotope 207 Bi. 17 18 12 A prior report existed finding 
that the amount of contamination largely depends on where the Bi ore comes from. 18 The BGO scintillators 
now commercially available have substantially reduced 207 Bi contamination.. 19 

Radioactive contamination (natural and induced) in the detection part contributes to the background 
more directly and should be absolutely minimized. Radioactive contamination has been measured for two 
high-light-yield scintillation materials with fast decay times, GSO(Ce) 5 6 and YAP(Ce) 20 21 (see Table 1 ), 
and their activation characteristics have been studied by irradiating them with protons of kinetic energies 
typical to the satellite orbit (~100 MeV). 10 12 22 23 The study has shown that the number of long-life line 
7-rays in the energy range of HXD is comparable for 2 cm thick YAP and 5 mm GSO for a fixed dose of 
proton. Note that the radiation length is quite different for the two crystals (see Table 1). We found recently 
that scintillation light yield decreases as temperature drops below 0°C for YAP but increases for GSO as 
shown in Fig. 5. 7 We also found that a slowly decaying scintillation component becomes dominant below 
0°C for YAP. 7 We plan to set the operating temperature at around — 15 ~ — 20° C so that the leakage current 
of the PIN diodes be reduced: hence GSO(Ce) became our final choice. 

We have already received about 40 GSO scintillators from Hitachi Chemical: their light yield are around 
30 % of a typical Nal(Tl) and give about A E(fwhm) ~ 6.8 — 7 . 5 % at 662 keV when directly, coupled to the 
PMT. When viewed through the BGO well, the resolution deteriorates typically to ~ 9 . 5 % at 662 keV at 
room temperature. Since the light yield is expected to increase by about ^ 15 % in the operating te mperatur e 
(— 20 ° C) as seen in Fig. 5, we expect the resolution to be around A E(fwhm) — 8.8% or ~ 7 %/ y/ EfMeVj. 


2.1.1 Fine collimator 

The two other experiments in the Astro E mission focuss on fine spectroscopy with imaging capability 
and have a narrow field of view (FOV) of around 17 x 17 arcmin 2 . The FOV of HXD is determined by 
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Figure 5: Temperature dependence of light yield for GSO, YSO and YAP when measured with a shaping 
amplifier with 500 ns time constant. 


Table 1: Characteristics of Nal, BGO, GSO and YAP 


Chemical composition 

Nal(TY) 
Nal (T£) 

BGO 

Bi4Ge 3 Ox2 

GSO(Ce) 
Gd 2 Si0 5 (Ce) 

YAP(Ce) 
YA10 3 (Ce) 

EfF. atomic number 

50 

74 

59 

35 

Density (g/cm 3 ) 

3.7 

7.1 

6.7 

5.5 

Rad. length (cm) 

2.6 

1.2 

1.4 

2.6 

Index of refraction 

1.85 

2.15 

1.9 

1.94 

At around 20° C 

Decay time (ns) 

-230 

—300 

—60 

-30 

Light yield (relative) 

100 

-12 

-28 

.—35 

Peak emission (nm) 

410 

480 

430 

347 

At around — 20° C: data on GSO and YAP are preliminary. 

Decay time (ns) 

-500 

—600 

—80 

—30 and > 500 

Light yield (relative) 

-75 

-15 

-30 

—20 

Peak emission (nm) 







Energy (keV) 


Figure 6: Limit to the sensitivity set by 
the diffuse background (B.G.L.) and source 
confusion (C. L.), with the 8x8 fine colli- 
mator (FOV=0.5) and without (FOV=4). 
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Figure 7: Transmission function of the fine collimator. 
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Figure 8 : Reflectivity of reflective sheets as 
a function of wavelength . 



Figure 9: Effective area of one of the 4 anit-counter 
planes used as the 7 -ray burst detector (BGO 2.0 cm) 
compared with that of BATSE one unit. 


the depth of the BGO well and the cell size of the fine colimator, both of which are in principle adjustable. 
In practice, however, the mechanical stability and availability of the BGO crystal set a limit at around the 
present value of 32 cm. The transverse dimension of the detection part is required to be substantially larger 
than the depth, setting the FOV at around 4.3 deg. At the low background rate of the well-type phoswich 
counter, the sensitivity will then be limited by the source confusion and diffuse background for the 4.3 deg 
FOV up to about 100 keV (see Fig. 6). 24 The passive fine collimator has been designed to reduce the source 
confusion and diffuse background . 24 

The collimator is made of phosphorus blonze sheets spot- welded to form 30 cm long 8 x 8 cells. The 
thickness of the phosphorus blonze sheets is 50/xm which gives a narrow FOV below 100 keV but secures an 
order of magnitude wider FOV for the 511 keV line (see Fig. 7). 

The challenge we now face is accuracy of the cell-to-cell alignment as well as the mutual alignment of the 
4 collimators in a detector unit. Possible nuclear activation of phosphorus blonze in orbit also concerns us 
and will be studied by irradiating the material with 100 MeV protons. 


2«1.2 UV-tansparent epoxy and reflector 


The shielding part (the well and the bottom in Fig. 2 ) is made of seven BGO plates and a BGO block: they 
are glued together by epoxy resin transparent in A ~ 350 — 600 nm, EpoTek 301—2 by Epoxy Technology 
Inc . 25 26 The detection part, GSO, is glued by silicon compound KE103 by Shin’Etsu Chemical 27 that 
can absorb the difference in the expansion coefficients of BGO and GSO (~ 10 -5 in one axis). It is also 
transparent in the wavelength range . 25 

The scintillator assembly is then wrapped by a light reflector sheet that has a high reflectance in the 
wave band and chemically stable. Our test has shown that the white GorTex sheet 28 has the. best reflectance 
as shown in Fig. 8 . 


2.2 Anti-Counters 

The anti-counters serve primarily to guard the 16 well-type phoswich counters from nuclear particle 
bombardment and to reduce the background in the phoswich counters (see Fig. 1 ). They are made of BGO 
and have wedge shapes as shown in the figure. Combination of these 20 anti-counters and the shielding parts 
shields all GSO scintillators with at least ~5 cm of BGO all around. This reduced the cosmic proton flux 
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on the GSOs by an order of maginitude . 23 The anti-counters also serve to reduce Compton scattered events 
as well as nuclear activation background events. 

The total area covered by these thick BGO anti counters ( 2.6 cm in average) is quite impressive: each 
face has a geometrical area of ~1200 cm 2 , and an effective area at 1 MeV of ^600 cm 2 as shown in Fig. 9. 
We plan to use all 4 faces to monitor 7 -ray bursts and transient phenomena. By comparing the counts in 
the 4 faces, we will be able to determine the direction to ~ 5°. 

Because of the scintillator shapes, we do not expect a good energy resolution in this monitor. The energy 
range coverage, on the other hand, will be quite wide as shown in Fig. 9. The scintillator used (BGO) is 
relatively fast and allows us to record transient phenomena in fine time bins (1/64 ms) as will be discussed 
in Sec.4,1. We believe this will be quite important when locating 7 -ray bursts by the arrival times registered 
by 3 or more satellites, one being Asiro-E. 


2.3 Silicon PIN Diode 

In the well-type phoswich counter, two layers of 2 mm thick FIN diodes sit in front of the GSO scintillators. 
Softer X-rays will be photo- asborbed in the two layers of PIN diodes, while harder photons pass through the 
two diode layers and reach the GSO crystal (see Fig. 3). 

The PIN diodes are introduced to fill the possible gap in the energy coverage between the phoswich 
counters (> 50 keV) and the CCDs (< 10 keV). The BGO wells provide the PIN detectors with a very 
low background environment and the diodes, in return, act as anti- coincidence shields for the scintillators 
against low-energy charged particles. The background level is expected to be substantially lower than 
10 - 5 /sec/cm 2 /keV for the PIN diodes. 

There are several technically critical issues in developing thick PIN diodes. One needs ultra high purity 
(~ ultra high resistivity) silicon wafers that give little volume leakage current . 8 29 30 Even with such wafers, 
extreme care must be taken in the diode fabrication process not to increase edge leakage current nor to lower 
the breakdown voltage . 31 32 Sample diodes with thickness of 1-1.5 mm have been produced by Hamamatsu, 
Micron, and Seiko Instruments Inc . 31 32 We plan to develope 2 mm thick diodes and stack two of them to 
obtain 4 mm effective thickness. To reduce the leakage current to a reasonable level (< a few nA), we plan 
to operate HXD at around — 20 ° C. 


2.4 Low Noise Charge Amplifier 

Improving the energy resolution is the most demanding issue in developing our silicon PIN diode. The 
energy resolution can be improved by optimizing the design of the preamplifier and the shaper amplifier. In 
reality, however, the total input capacitance (C in ) which is the sum of the junction capacitances of two PIN 
diodes (~ 2 x 20 pF) and the capacitance of the cables connecting the diodes to the preamplifiers (~ 30 pF) 
will determine the resolution . 8 29 30 We currently set our goal at an energy resolution of AE(fwhm) ~ 3 
keV at — 20 ° C and are developing a low-power low-noise amplifier system optimized to this large input 
capacitance. The current design gives the equivalent noise figure shown in Fig. 10 . 


3 HOUSING OF THE DETECTOR ARRAY 

The mechanical vibration and shock is expected to be much higher for the M— V rocket of ISAS than for 
rackets powered by liquid fuel and the housing for the detector array requires great care. The dissipation 
of the heat generated by the AE parts placed in the housing (~ 15 W) poses additional constraint to the 
design. Since the BGO crystal is quite brittle, some support mechanism must be installed to damp external 
vibrations and to minimize the internal stress of BGO. The thermal expansion coefficient of BGO is an order 
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Figure 10: Equivalent noise count (ENC) of Figure 11; Biock diagram G f the HXD electronics 

the preamplifier for PIN diodes. 


of magnitude smaller than any metal and consequently some critical parts will be made of carbon fiber or 
glass fiber reinforced plastic composite. 

The detector array will probably be serviced including replacement of some detector units. This requires 
each unit to be mechanically independent. On top of all the above, the total amount of passive material 
must be minimized because of background generated by bombardment of cosmic protons and neutrons. 


4 ELECTRONICS AND DATA ACQUISITION 

The electronic system for the Astro-E HXD consists of the analog electronics part (AE) and the digital 
electronics part (DE) 9 as shown in Fig. 11. The AE consists primarily of pulse shape discriminators (PSD) 
and circuits associated with ADCs. The DE formats and processed the digitized data. The final data are 
then transmitted to the ground stations under control of DE. Data from various environmental monitors will 
also become critical in the HXD analysis, because the scintillation light yield is known to be temperature 
dependent and reliable background subtraction requires accurate monitoring of possible background origins. 


4.1 Analog Electronics 

The AE consists of one control board (ACU) and eight signal-processing boards (see Fig. 11). The latter 
are broken into the following two types: WPU for the well-type phoswich counters and PIN diodes, and 
TPU for the anti counters. The ACU handles power supply to AE and monitors the house keeping data 
about power supplies, timing, temperatures, and others. 

The analog signal is taken from the Iasi or the second last dynode of the PMT in the well-type phoswich 
counter. It is first fed to a charge sensitive preamplifier placed in the detector housing and then to a pulse 
shape discriminator (PSD) in the WPU board. The anode signal of the PMT is used for a fast pre-trigger 
that generates the peak-hold gate for the PSD. The PSD selects signals whose time profiles are consistent 
with that of the detection part (GSO), and rejects those contaminated with slow-decaying BGO scintillation 
light. While in orbit, the counting rate in the shielding part will be 4-5 orders of magnitude higher than 
that of clean hits in the detection part. The PSD selects the clean hits by comparing the output pulse 
heights from the two shaping amplifiers with different integration times. The PSD has been developed as a 
semi-customed LSI to reduce the circuit size and power consumption. 33 34 35 

PIN diodes can also trigger the data aquisition system if a diode level discriminator detects a signal while 
no slow scintillation light is detected by the PSD for the phoswich unit housing the diode. 
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Figure 12 : Correlation between the fast and slow 
shaper outputs (r = 100 ns and 500 ns) obtained 
by the EM unit at — 20 ° deg. and the PSD LSI 
irradiated by 662 keV gamma-rays. See text. 



Figure 13: Pulse height distribution of the slow 
shaper output selected by the PSD gate and of- 
fline filtering from Fig. 12. Since the data was 
taken with only one unit and in a temperature 
control box, gamma-rays escaped through the 
3 mm BGO wall after hitting GSO and stray 
gamma-rays entered through the 3 mm BGO wall 
to hit GSO. 


When a trigger is issued either by a PSD or a PIN diode, all analog outputs from the corresponding 
phoswich unit (those of the slow shaper, the fast shaper, and the four diode shapers) are digitized. The 
upper discriminator outputs (see below) and existence of any hits (the hit pattern) of all the 16 counters 
and 20 anti-counters are recorded together with the time of the event by the WPU board. 

When large inorganic scintillators are used in the space as in the present case, highly ionizing nuclear 
particles often deposit energy higher than the saturation point of the amplifying system. They may poten- 
tially cause some erroneous response in the analog system. In the AE, we install discriminators to detect 
such events (the upper discriminators) and protect the analog system from such misoperation. 

The signal from the 20 anti-counters are processed in units of five in the TPU boards. The outputs 
of all 20 lower discriminators are distributed to the WPU boards via the “hit pattern” bus. Signals are 
then summed to produce two sets of pulse height distributions: one has only 4 energy bins but renewed 
every 1/64 msec, and the other has 64 energy bins and renewed at 512 msec interval. These transient pulse 
height distributions are used to detect 7 -ray bursts, various flare-ups, and other bright transient phenomena. 
Once a 7 -ray burst occurs, for example, the on-board burst detection circuit or the program running in the 
CPU triggers the memory control circuit of the histogram buffer to record the time-resolved pulse height 
distributions of the burst for 128 seconds. We expect to detect about 160 bursts per year. 


4.2 Pulse Shape Discriminator 

The PSD distinguishes events from the two scintillators, the fast-decaying scintillator (GSO) and the 
slow-decaying scintillator (BGO). There are several PSD methods developed and used in 7 - ray astronomy 
and nuclear physics, one of which being the double shaping method adopted in the HXD. In this method, 
the signal from a phototube is integrated with two different shaping times and the two output pulse heights 
are compared. This method is less affected by external electronic noise, while its demerit is that the cir- 
cuit becomes somewhat complicated. This demerit can be solved by implementing the circuit on an LSI 
chip . 33 34 35 In our LSI chip, the two integration times are set at around Tf ast = 100 ns and r s [ ow = 500 ns. 
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The outputs of the two shapers have different pulse heights for BGO signals, while they are almost equal 
for pure GSO signals. The pulse-shape discimination LSI has been developed using a semi-customed LSI 
technology . 33 34 35 Shown in Figs. 12 and 13 is the result of a test done by using the LSI circuit with the 
phoswich unit of the Engineering Model: the irradiated gamma-rays are from 137 Cs (662 keV). In Fig. 12 , 
the linear cluster extending at an angle around 30 deg. to the ordinate, is due to hits on the BGO shielding 
part while the one extending to right is due to clean hits on the GSO detection part. Compton scattered 
gamma-rays form two linear clusters bridging between the blob corresponding to the clean hits on GSO and 
either one of the two blobs corresponding to the bottom BGO and the collimator BGO. One can see that 
pure GSO events can be seperated from BGO events and Compton scatterd events. Fig. 13 shows how the 
background is reduced by selection of the GSO cluster in Fig. 12 by the PSD chip and subsequent off-line 
analysis. Since we have put only one unit in a temperature- controlled box, Compton- scattered gamma-rays 
escaped or entered through the 3 mm BGO wall as apparent in Fig. 13. In the real detector, these 7 -rays 
will be captured by neighboring units and rejected when the hit pattern is interrogated. 


4.3 Digital Electronics 

The DE controls acquisition of the data, coordinates among them, detects and reacts on requests for 
services, and provides the primary interface with the satellite data processor for command and telemetry 
(see Fig. 11 ). These tasks are executed by a system consisting of one Intel 80386 running at 8 MHz. To 
obtain a high data- acquisition, rate, each AE board is capable of sending data by the DMA transfer mode 
to the DE. Here the data for 64 events are blocked so that the interrupt rate to CPU remains low . 9 

The CPU also performs some software event selections by examining the hit pattern around the unit in 
which the trigger is initiated: this further reduces Compton scattered events and particle interaction events. 
The CPU acquire calibration data and data used to estimate background. 


5 CALIBRATION AND MONITORING 


On board calibration will be needed in the energy measurement by the PIN diodes and phoswich detectors. 
For the diodes, we are examining possiblity of using K X-rays from the phosphorus blonze of the fine 
collimator. For GSO, the o-rays emitted by a natural isotope 152 Gd, appearing at around E 7 ^ 385 keV 
with a counting rate around ~ 10 ~ 4 /s/cm 2 /keV, will be used . 10 11 12 13 14 

Critical items to be monitored include various counting rates that will be needed in estimating the 
activation background. The temperature of the scintillators and phototubes are also to be monitored because 
of their temperature dependence shown in Fig. 5. 


6 EXPECTED BACKGROUNDS 


Possible in-orbit backgrounds have been estimated for the present counter design: GSO of 2.4 x 2.4 x 0.5 
cm 3 shielded by BGO of ~ 5 cm thickness all around . 22 8 The proton flux has been assumed to be that of the 
model given in the reference 23 at the solar minimum, although the scheduled flight time corresponds roughly 
to the solar maximum when the background is expected to be about half . 23 The results are summarized in 
Table 2 : one can see that the internal and cosmic-ray induced radioactivities will be the dominant source 
of background below 300 keV. We therefore anticipate that the sensitivity of HXD will be limited by these 
induced radioactivities. 

Among the listed background origins, prompt backgrounds produced by charged particles on BGO can 
easily be rejected. Gamma-rays generated in passive material around the detector array have some chance 
to leak through a few centimeter of BGO, if their energy is above ~ 300 keV. We have to minimize passive 
material close to the HXD detectors, eg. the HXD housing and the fine collimator. Gamma-rays generated 
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Table 2 : Expected background rates 


In-orbit activation 
Leak-thru 7 -rays 

Off- aperture CXB 
Radioactive impurities 


~ (0.5 ~ 1) x 1G -5 /sec/cm 2 /keV at 100 keV 
~ 1 x IQ -6 /sec/cm 2 /keV at 300 keV 
~ 2 x 10~ 5 /sec/cm 2 /keV at 400 keV 
~2 x 10 -4 /sec/cm 2 /keV at 500 keV 
~ 7 x 10 -6 /sec/cm 2 /keV at 50 keV 
~ 1 x 1 G -6 /sec/cm 2 /keV at 100 keV 
5^1 x IQ -6 /sec/cm 2 /keV below 300 keV 


by nuclear interactions in the tip of the XRT complex have some chance to enter the FOV determined by 
the RGO well. 

What may become a dominant background and yet is difficult to estimate is faking of clean-hit GSO 
events by pile-up of BGO hits or by Cherenkov photons in the PMT face windows. Background due to 
radioactivity and amount of radioactivity that may be produced have been and still are actively studied by 
bombarding the materials around the HXD detector with protons and by running the simulation programs. 22 
Through these efforts and adequate on board monitoring, we expect to be able to predict the background 
for actual observations to an accuracy better than ~ 5%. We note that such predictions worked to about 
~ 1% level in the Ginga-LAC. 16 


7 CONCLUSIONS 


The design and expected performance of the Hard X-ray Detector (HXD). one of the 3 experiments 
aboard Astro-E have been presented. Most of its critical elements have been prototyped, tested, and verified 
their validity. One item that still requires continuing R/D works is the 2 mm thick silicon FIN diodes. The 
total weight and the total power consumption are to be reduced further. We expect the above issues to be 
solved within this calendar year. 
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